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The major and minor axes of the polarization elhpses that surround singular lines of cir- 
cular polarization in three dimensional optical ellipse fields are shown to be organized into 
Mobius strips. These strips can have either one or three half- twists, and can be either right- 
or left-handed. The normals to the surrounding ellipses generate cone-like structures. Two 
special projections, one new geometrical, and seven new topological indices are developed to 
characterize the rather complex structures of the Mobius strips and cones. These eight in- 
dices, together with the two well-known indices used until now to characterize singular lines of 
circular polarization, could, if independent, generate 16, 384 geometrically and topologically 
distinct lines. Geometric constraints and 13 selection rules are discussed that reduce the 
number of lines to 2, 104, some 1, 150 of which have been observed in practice; this number 
of different C lines is ~ 350 times greater than the three types of lines recognized previously. 
Statistical probabilities are presented for the most important index combinations in random 
fields. It is argued that it is presently feasible to perform experimental measurements of the 
Mobius strips and cones described here theoretically. 



I. INTRODUCTION 



Three dimensional (3D) optical fields are, except in special cases, elliptically polarized. In 
paraxial fields the polarization ellipses lie in parallel planes oriented normal to the propagation 
direction, but in 3D fields the ellipses generally have a wide range of spatial orientations. 

The generic singularities of 3D ellipse fields are meandering lines of circular polarization, C 
lines, and meandering lines of linear polarization, L lines [1 — 16]. Here we study the arrangement 
of the polarization ellipses surrounding C lines; the arrangements of the polarization ellipses that 
surround L lines will be reported on separately. Where a C line pierces a plane, S, a point of 
circular polarization, a C point, appears. Previous studies of C lines have concentrated on the 
two dimensional projections onto S of those ellipses whose centers lie in this plane on a circle that 
encloses the C point. Examining the full 3D arrangement of these ellipses, we find that their 
major and minor axes generate Mobius strips with either one or three half- twists; examples of such 
strips are shown in Fig. [TJ 




FIG. 1: Mobius strips surrounding C lines in a computed random optical field. The strips shown here are 
generated by the major axes of the ellipses located on a small circle that encircles the line. The Mobius 
strips in (ai)— (as) contain a single half- twist, those in (bi) — (bs) three half- twists, (ai) and (bi) are overall 
views of the strips that show their 3D structure. Here ellipse centers are shown by small gray dots, the 
major axes of the ellipses by thin straight lines, and the endpoints of these axes by small black dots. {0,2) 
and (b2) show the circular strips in (ai) and (bi) cut open and straightened to better display the twist 
structure. Here the ellipse centers form the thick central line, and, for clarity, only one half of the ellipse 
axes are shown (thin straight lines). The small white (black) circles mark the endpoints of the first (last) 
ellipse, and serve to emphasize the twist structure; when the strips are closed these ellipses are adjacent, 
but their endpoints lie on opposite sides of the circle of ellipse centers, (aa) and (bs) are views from above 
that illustrate the fact that these strips, which have an odd number of half-twists, have not only one side, 
but also have only one edge. Here the centers of the ellipses are shown by gray circles, the endpoints of 
ellipse axes that lie above (below) the plane of ellipse centers by white (black) circles. Starting at the 
highest point on the "top" edge (white arrow) and traversing a 360° circuit along this edge one arrives at 
the "bottom" edge (black arrow) . Continuing in the same direction along a second 360° circuit one returns 
to one's starting point on the "top" edge, showing that there is, in fact, only a single edge to the strip. The 
ellipse endpoints in (ba) can be seen to form a trefoil knot, whereas those in (as) do not form a knot. 



The three orthogonal principal axes of the (always planar, [17, Sect. 1.4]) polarization ellipse 
are its major and minor axes, and the normal to the ellipse. On a C line the major and minor 
axes of the ellipse become equal, and the ellipse degenerates into a circle, the C circle. Because a 
circle has no preferred direction, its major and minor axes are undefined (singular). The normal 
to the C circle, however, remains well defined, as do all three principal axes of the ellipses that 
surround the C line. The projections of the major and minor axes of the ellipses whose centers 
lie in a plane S pierced by a C line rotate about the C point, generically with winding number 
(net rotation or winding angle divided by 2tt) I = ±1/2 [1,3, 10]. Examples of such rotations are 
shown in Fig. [2] for the major axis; in all known cases, for a given C point the winding number of 
the projections of the minor axes onto S is the same as that of the major axes [1, 3, 10] 
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FIG. 2: C point winding number /. Shown in panels (a) and (b) ((c) and (d)) are C point Mobius strips with 
one (three) half-twists, upper panels, and their projections onto the xy-plane, lower panels. As in Figs, la 
and lb, in the upper panels ellipse centers are shown by small gray dots, major axes of the ellipses are shown 
by thin straight lines, and the endpoints of these axes by small black dots. The central black dot is the C 
point. In the lower panels the ellipse centers form the black circle, and the projections of the ellipse axes 
are shown by short straight lines. During one 360° circuit around the C point in say the counterclockwise 
direction, in (a) and (c) ((b) and (d)) these lines rotate by 180° in the same counterclockwise (opposite, 
clockwise) direction as the circuit, and / = +1/2 (—1/2). 

Winding number / and a geometrical index described later, the line index [1, 3, 18], fully char- 
acterize the projections in Fig. [21 together these indices lead to the currently known three distinct 
C lines [1,3]. As is evident from Fig. [21 however, these projections of the major/minor axes are 
insufficient to determine the properties, or even the existence, of the parent Mobius strips. In 
later sections we introduce one new geometrical, and seven new topological indices to characterize 
these strips; together these eight new indices increase the number of distinct C lines by a factor of 
~ 350. 

The plan of this report is as follows: In Section II we describe the analytical and numerical 
tools used in the computer simulations employed in later sections to study 3D ellipse fields, their 
C lines, their Mobius strips, and related structures. These tools are similar to those we used to 
study the one-full-twist Mobius strips that surround ordinary (i.e. nonsingular) ellipses [14, 15]. In 
Section III we describe in detail the Mobius strips and related structures that surround C lines, 
introduce the new indices that characterize these lines, and extend the important line classification 
of C points [1,3, 18]; although this classification does not involve an invariant topological index, it 
does serve to further characterize the arrangement of the ellipses that surround these points. In 
Section IV we present statistical data for the relative occurrences in random fields of the many 
different C lines. We summarize our main findings in the concluding Section V. 



Coherent measurement techniques permit the determination (amplitude and phase) of all three 
orthogonally polarized components of 3D microwave fields [6,7,19,20]; recent advances in inter- 
ferometric nanoprobes provide similar capabilities for optical fields [21 — 28]. It is therefore now 
possible to carry out experiments that can measure the highly unusual structures described here. 
C lines are degeneracies of the matrix that describes the 3D polarization ellipse [9, 29, 30]. Other 
physical systems such as liquid crystals, strain fields, flow fields, etc., are described by similar matri- 
ces whose degeneracies can be expected to yield analogs of C lines that are likely to be surrounded 
by Mobius strips. 

At present there are no known practical applications for optical Mobius strips; as these strange, 
engrossing objects become better understood, however, useful applications may emerge. 

II. METHODS 

As indicated in the Introduction, the only well defined direction associated with a C point is 
the normal to the C circle. We call the plane that is perpendicular to this preferred direction, i.e. 

the plane of the C circle, the principal plane, and denote this plane by So- We emphasize that 
So need not be, and in general is not, perpendicular to the C line, and that the orientation of the 
principal plane relative to the C line changes as one moves along the line [1, 3, 8 — 10, 13 — 15]. 

A. Ellipse Axes 

Here we describe the methods we use to calculate the axes of the polarization ellipses that 
generate the Mobius strips surrounding C points and C lines. 

We label the major and minor axes of, and the normal to, the general polarization ellipse by ot, 
f3, and 7, respectively. Given an expression for the (here complex) optical field E, a, /3, and 7, 
can be calculated in two seemingly different ways. The first, due to Berry [8, 10], is 

a = Re(EVE • E), (la) 
p = Im(EVE-E), (lb) 
7 = Im(E*xE). (Ic) 

The second involves finding the three eigenvalues, Aj, and three normalized eigenvectors, i/j, z = 
1, 2, 3, of the 3x3 real coherency matrix, 

Mij = Re{E*Ejy, i,j = X, y, z. (2) 

The largest eigenvalue, here A3, is associated with the major axis of the polarization ellipse, the 
next largest eigenvalue, A2, with the minor axis, and the smallest eigenvalue Ai = 0, with the 
normal to the ellipse. 

Ai = because for the monochromatic fields assumed here the polarization ellipse is always 
planar [17]. In terms of A, the length a of the major axis of the polarization ellipse is given by 
a = -v/As, the length b of the minor axis by 6 = s/)^. 

The above two seemingly different methods are reconciled by noting that a/ |q| = 1/3, /3/ |/3| = 
1/2, and 7/ I7I = t^i- Although an analytical proof of these equalities is still lacking, we note that 
in every one of the literally hundreds of cases studied these equalities were found to hold to within 
numerical accuracy. In what follows, a, /3, and 7 are, without change of notation, understood to 
each be normalized to unit length. 

As we move through the wavefield, in order to ensure that a and f3 are smooth, continuous, 
single valued functions of position, we calculate ^E^ as follows: We write E • E =E^ + E^+E'^ = 



Aexp (iip), unfold (unwrap) if as needed to eliminate spurious discontinuities of 27r, and then write 
VE • E=\/]4exp {i(p/2). 

Although the polarization ellipse has inversion symmetry, 7 in Eq. px|) is an axial vector that 
defines a unique positive direction. This direction, which is determined via a right hand rule from 
the way E traces out the polarization ellipse in time, permits us to define a unique orthogonal xyz 
coordinate system in which the positive z-axis is along the positive direction of 7, the x and y axes 
are along a and /3, and xyz form a right-handed 3-frame. 

C lines may be located in two different ways: The first is to locate the zeros of E • E [8 — 10]. 
The second is to locate the zeros of the C point discriminant Dc = af — 4a2', this discriminant is 
obtained from the characteristic equation + aiA + 02 = of the real coherency matrix Mij in 
Eq. ([2]) [29]. As expected, both methods were found to yield the same result. 

B. Axis Projections 

Having traced out a C line, we move along the line to some C point and find 7 for that point. 
As noted above, 7 for a given point is the positive normal to Sq at the point. In looking at the 
ellipses in Sq we will always look down the positive z-axis (i.e. +z points towards us). 

In constructing the Mobius strips and their projections we focus on those ellipses whose centers 
lie in T,q on a small circle that surrounds the C point. We call these ellipses the surrounding ellipses 
and label the corresponding circle (Tq. For convenience, we take the center of ctq to coincide with 
the C point, although any other simple path in Eq centered on the C point yields the same result 
for the various topological indices as does (Tq. 

In addition to projecting the ellipse axes onto Sq, we study two other projections. The first, 
and most important, we call the tq projection. In this projection we erect a rotating plane, the tq 
plane, tq is oriented normal to Sq, contains the normal to the C circle, axis 7 (the +z axis), and 
rotates around this axis. Where do pierces the plane we establish an orthogonal XoZg-coordinate 
system. The angle of rotation xo of the plane tq is measured counterclockwise from the fixed 
(laboratory) x-axis in Eq. Projected onto tq is the axis of the ellipse whose center lies in ctq at 
the point where tq intercepts this circle; we label this point Pq. Pq, of course, moves along ctq as 
Xo increases — 27r. Fig. ^ displays the relevant geometry. 

The second projection we call the ttq projection. In this projection we erect another rotating 
plane, the ttq plane. Like tq, ttq is oriented normal to Sq, but unlike tq, ttq is tangent to fio. We 
establish an orthogonal YoZo-coordinate system at the point of tangentsy, the point Pq, such that 
the triplet XqYqZq forms a right handed coordinate system. The same angle, xo; that measures 
the rotation of tq measures the rotation of ttq. Projected onto ttq is the axis of the ellipse whose 
center lies at Pq. Fig. [8)3 illustrates the geometry of this case. 

Thus, in determining the geometrical and topological properties of the 3D Mobius strips formed 
by a given axis we use three orthogonal projections — the minimum number of projections required 
for the characterization of a three-dimensional object. 

In general, there are therefore nine different projections: the projections onto the three planes 
S, T, and vr, of each of the three axes a, (3, and 7 of the ellipses on a. Although these projections 
are interconnected, as becomes apparent, they yield a multitude of different topological winding 
numbers This may appear surprising because, after all, the three-frame afSy has only three 
rotational degrees of freedom. But the properties of the Mobius strips depend on the rotations of 
all the ellipses on (Tq. These ellipses are unbounded in number, and their three- frames can rotate 
independently, subject only to the restrictions of continuity. A practical, albeit not fundamental, 
further restriction is introduced by the fact that sufficiently close to the C point a linear expansion of 
the field suffices to determine what happens in the generic case. As becomes apparent, this leads 



to additional interconnections between indices, and forbidden combinations of indices (selection 
rules), so that structures that are allowed in principle by geometrical and topological constraints 
may not appear in practice. 
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FIG. 3: Projections tq and ttq. Shown as an example is axis a projected onto plane tq (plane ttq) in (a) 
in ((b)). Other symbols are explained in the text in Section IIB. Position on ctq is parameterized by the 
angle xoj or by arclength sq = roXo, where ro is the radius of ctq. 



C. Computed Optical Fields 

1. Random Speckle Field 

We study here two types of computer simulated random optical fields. The first is a field 
composed of a large number of randomly phased, linearly polarized plane waves with random 
polarization and random propagation directions. This speckle field (speckle pattern), described in 
full detail in [13 — 15], has the advantage that it is an is an exact solution of Maxwell's equations. 
It has two disadvantages, however: The first is that C lines must be found and followed empirically, 
using a time consuming search based on either the zeros of E • E, or the discriminant Dc of the 
coherency matrix described above. Thus, the number of different lines that can be studied is 
limited, and as a result, a limited menu of different Mobius strips is available for detailed study. 
Nonetheless, this random field serves as our gold standard, and is used to verify the existence of 
the structures measured with the aid of the linear expansion described below. 

2. Linear Expansion 

Locating the origin at an arbitrary point in a generic complex optical field, the three Cartesian 
components of the field in the immediate vicinity of the point can be expanded to first order (the 
order required here) as 

■^x \ I Ojx "I" i^a; \ / Pxx "I" ^Qxx Pxy "I" '^Qxy Pxz ~^ ^Qxz \ f \ 

Ey = ay+\hy + Py.j; + \Q Pyy + iQ yy Py ^ + \Q y Z ^ ' 

Ez ) \ O-z ~l~ '^z j \ Pzx ^ ^Qzx ^zy + '^Qzy Pzz + / \ ^ / 

where, for example, Pxy = ^e{dEx/ dy)^, Qxy = lm{dEx/dy)Q, etc. The divergence condition 
V • E = can be satisfied by setting Pzz = — {Pxx + Pyy), and Qzz = — {Qxx + Qyy)- 



At the selected point, which can have arbitrary ehiptical polarization, we erect a convenient 
coordinate system in which the +z-axis is along the positive normal 7 to the plane of the ellipse, 
the X-axis is along the major ellipse axis a, and the y-axis is along the minor axis (3. In this 
coordinate system the plane z = is the principal plane Sq of the ellipse at the origin. 

When the point at the origin is a C point on a C line, the condition (E • E)q = can be satisfied 
by setting ay = = 0, = bz = 0, and = it \ax\- Thus, in Sq all constraints are satisfied 
without loss of generality by writing for a C point 

Ex = a + {Pxx + '^Qxx) X + {Pxy + iQxy) y, (4a) 
Ey = is \a\ + {Pyx + iQyx) x + {Pyy + iQ yy) v, (4b) 

Ez = {Pzx + '-^Qzx) X + {Pzy + iQzy) V, (4c) 

where s = ±1, and a and the P and Q remaining in Eq. (jl]) are unconstrained. 

Simpler expansions are also possible. The one-half-twist Mobius strip in Fig. [1^ is closely 
approximated by Ex = — 1 — (1 — i) y, Ey = —i + x, Ez = —ix + y, the three-half-twist strip in Fig. 

hy Ex = I - y,Ey = i - x,Ez = y. 

We make contact with the random speckle field described above by generating the various 
constants in Eq. (jlj) with the aid of the procedure described below: this procedure is appropriate 
to a random field whose components obey circular Gaussian statistics, and ensures that the Mobius 
strips we study have generic properties that are likely to correspond to those in real physical fields. 

We start by writing the joint probability density function (PDF) *p of the three Cartesian field 
components in Eq. ([3]), which describes a general point in the fixed laboratory xyz frame, as 

^ {Ex, Ey, Ez) = ^ {Ex) ^ {Ey) <p {Ez) , (5a) 

V{Ej) =V{aj,bj,Pjx,:;Qjz) , j = x,y,z (5b) 

*P {aj, bj,Pjx, Qjz) = Vg {aj) *Pg {bj) {Pjx) -Vg {Qjz) , (5c) 

1 / u'^ 



But we want {Ex, Ey, Ez) not at a general point, but at a special point, a C point, and not in 
the fixed laboratory xyz frame, but in the principal axis x'y'z' frame of the C point. Accordingly, 
we proceed numerically as follows. In accord with Eq. ([5]) we first choose the sign of s and 
the values of the various a, b, P, and Q in Eq. ([3]) by consulting a random number generator 
that produces a Gaussian distribution with, for convenience, unit variance. We then adjust the 
parameters a and b so as to satisfy (E • E)q = 0, writing and bz in terms of ax, ay, bx, by. Next, 
we find the eigenvectors of the coherency matrix in Eq. ([2]); these eigenvectors yield the direction 
cosines of the principal axes x'y'z' of the C point relative to the laboratory xyz frame. As discussed 
above, eigenvector vi is the well defined normal to the C circle (axis f), whereas i'2 and 1^3 are 
arbitrary orthogonal directions in the plane of the C circle. 

We form a matrix D from these direction cosines, and use it to transform the parameters a, 
b, P, and Q, in the laboratory frame to a parameter set a', b', P' , and Q' in the principal axis 
frame of the C point. Forming a vector a (vector b) from the a (6) parameters, and matrices P 
and Q from the P and Q parameters, respectively, and transforming both coordinates and field 
components, we obtain 

a' = D • a, (6a) 
b' = D • b, (6b) 
P' = D P D \ (6c) 
Q' = D Q D 1. (6d) 



As expected, we find a'^ = b'^ = 0. Writing a',j.+ \b'^ = exp (i(/?^), a'y+ib'y = ^4^ exp (i^?^) , 
we find, again as expected, A'^ = Ay and \ip'^ — ipy\ = 7r/2. But this is not yet in the convenient 
form of Eq. We complete the transformation by writing E" = exp(— i6^/a^)E', extract the 

real and imaginary parts of E", set z = 0, and after dropping for convenience the double prime 
superscripts obtain Eq. 

The individual PDFs of P' and Q' remain circular Gaussians with unit variance. This is 
not surprising, because P and Q are transformed under a similarity transformation in which the 
arguments of the cosines in D are uniformly distributed — 2tt. The above procedure is required 
because it preserves specific correlations between the final P" and Q" that are absent in the initial 
P and Q; these correlations arise because the elements of D are themselves functions of P and Q 
[15]. 

In contrast to the PDFs of P and Q, the PDF of a, shown in Fig. 4, differs markedly from a 
Gaussian. 
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FIG. 4: Final PDF of a in Eq. Q. 

Although Eq. Q appears to describe an isolated C point at the origin, this C point is, in fact, 
part of a continuous C line. In addition to the C point at the origin, Eq. ^ can also produce 
one or more parasitic C points (C lines) at locations that depend upon the various parameters. 
As long as the surrounding circle do is sufficiently small that it does not contain, or come close 
to, these parasitic C points they cause no trouble. Occasionally, with a probability of ~ 1 — 2%, 
a parasitic C point is so close to the C point at the origin that it strongly distorts the field in 
the region of the central point. When this happens the Mobius strip at the origin can become so 
distorted that its parameters cannot be accurately measured. In Section IV where we discuss the 
statistics of the strips we eliminate such cases (outliers) from the analysis. 

D. Graphics 

We discuss here three important aspects of the graphic representation of the Mobius strips and 
other wavefield structures that are presented throughout this report. 



1. Length scales 



We start by noting that there is no geometric or, indeed any, relationship between the radius 
of ctq, which is a true length, and the length of the line used to represent an ellipse axis. 

The length scale of (Tq is set by the variances of the Gaussian distributions in Eq. ([5]) that 
are used for the constants a and 6, which measure the field amplitude, and the derivatives P and 
Q, that measure the change in field strength per unit length. We take the field and derivative 
variances to be unity, which fixes the length scale of ctq to also be unity. When the radius rg 
of (jQ is sufficiently small (i.e. small compared to one), the field scales uniformly with tq and its 
structure becomes independent of this radius. Here, we take tq = 10~^, having verified that this 
is "sufficiently small" in the above sense. 

Although the length scale of the field is also unity, the length scale of the polarization ellipse 
remains arbitrary. The reason is that the ellipse is a representation of the magnitude and direction 
of the electric field vector over an optical cycle at a point. We can therefore choose any convenient 
scale to represent the ellipse, because in the physical field two arbitrarily close ellipses do not 
overlap. Thus, the widths of the ribbons representing the Mobius strips in Figs. [HandEl and in 
all other figures presented here, are adjusted as needed for clarity. Of course, in a given figure a 
single length scale is used for all ellipses. 

2. Scaling 

By continuity, the tilt of the plane of the ellipses on ctq relative to the plane Sq, the plane of 
the C circle, decreases to zero as the radius tq goes to zero. This tilt is of order rg, which for 
ro = 10"^ is of order 0.005". In later figures (such as Figs. [7] — I12|) we will show the projections 
of the ellipse axes onto the XZ-plane or the yZ-plane of Fig. [3l i.e. the planes tq and vrQ. In 
such figures we use anisotropic scaling, expanding the Z-axis relative to the X- or Y-axis so as 
to make visible the relative tilts of the axes of different ellipses, the absolute values of these tilts 
being unimportant in these figures. 

3. Paths 

As mentioned, we choose ao to be a circle centered on the C point for simplicity — other simple 
paths such as ellipses centered on the C point yield the same values for the topological indices as 
does fTQ. The path must be "centered " because when shrunk to zero it must enclose the C point 

The choice of a circle centered on the C point yields the most symmetrical possible form for 
the Mobius strips, enhancing our ability to understand their sometimes complex structures. But 
even within the domain of circles there is nothing special about a radius of say tq = 10~^ — 
ro = 1.1 X 10~^, or any other value, is obviously equally good. From this follows that the C point 
is surrounded by an infinite set of centered nested Mobius strips. 

It is, of course, difficult, if not impossible, to visualize the full 3D structure of such a field. 
Below we dissect out a single representative Mobius strip of the nested strips that surrounds C 
points, the do Mobius strip, and proceed to study this strip in detail. 

III. MOBIUS STRIP INDICES 

Throughout this report we restrict ourselves to the case in which the plane of observation S 
coincides with the principal plane Sq. In this plane we find that all indices are the same for axes 



a and /3, even though the detailed geometries of the Mobius strips generated by these axes are 
different. For simphcity and uniformity, in what follows all examples presented are for axis a. As 
S rotates away from Sq a complex set of phenomena set in: the universal equivalence between 
axes a and (3 is broken, right-handed Mobius strips transform into left-handed one, and vice versa, 
one-half-twist strips change to three-half twist strips, and vice versa, indices abruptly change sign, 
etc. These phenomena will be reported on separately. 

A. Topological Indices of the Projection onto Sq 

1. Indices Ia,p 

As indicated in the Introduction, the sole topological index used in previous studies to charac- 
terize C points and C lines is based on the rotation about the C point of the projections onto S 
of the major or minor axes of the ellipses on the surrounding circle a. When S = Sq (so that 
a = o"o) it is easily seen that index / (Fig. [2j) is the same for both axis o: and axis (3: sufficiently 
close to the C point the surrounding ellipses are tilted negligibly out of the plane Sq , the plane of 
the C circle, so that the projections of a and /3 onto Sq are very nearly orthogonal at all points on 
o"o. Thus, these projections are locked together, and as one rotates so does the other: as a result, 
la = If3- This is true even when T, ^ Sq. As expected, we find that in our simulations la = 1/3 in 
all cases (S = Sq, and S ^ Sq): accordingly in what follows we write this index as Ia,f3- 

Generically, la^fs = ±1/2 [1,3]. This is not due to geometric or topological constraints, but 
is rather a consequence of the linear expansion; when higher order terms dominate, higher order 
values for this index are possible [32]. Half integer values for /^^^ can occur because of the 
symmetry of the ellipse which returns to itself after rotation by vr; for the rotations of a vector, for 
example, only integer values for the index are possible. 

Hidden in the diagrams that illustrate la, (3 lie long overlooked clues to the fact that these 2D 
diagrams are projections of 3D Mobius strips. These clues are discussed in Fig. [5] 

2. Index /-y 

What about the projection onto T,q of axis -y of the ellipses on uo? As shown in Fig. El these 
projections rotate about the C point with integer winding number = ±1. 1^ = ±1 because 7 
is a vector. But wait! Axis 7 is well defined at the C point, and normally we associate winding 
numbers with singularities, i.e. with properties that are undefined; for example winding number 
la, 13 is associated with the fact that axes a and (3 are undefined at the C point. So is ly a valid 
topological index, or not? 

In answering this question in the affirmative we note the following: Although axis 7 is well 
defined, its projection onto Sq is singular. The reason is that although the projection of this axis 
for all the surrounding ellipses are lines with well defined directions, the projection of axis 7 for 
the C point itself is a point whose direction is undefined. Thus, in Sq, axis 7 and axes a and (3 
behave similarly, both projections have undefined directions — axes a and (3 because the C circle 
has no preferred direction, axis 7 because a point (a circle of vanishingly small radius) also has no 
preferred direction. 

The above comparison is closer still if instead of plotting the projections of the ellipses we plot 
the quantities given in Eq. ([1]). As defined there, a and f3 go to zero at a C point, and maps of 
the projections of a, f3, and 7, are similar — all three maps show a point at the location of the C 
point surrounded by lines that rotate about the point with winding number of ±1/2 (±1) for axes 
a and /3 (axis 7). 
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FIG. 5: Hidden clues. As shown in Fig. [21 in determining la.fi one traditionally focuses on the rotation of 
the projections onto Sq of the major or minor axes of the ellipses on the surrounding circle ctq- It is amusing 
to note in retrospect that these diagrams contain an important clue overlooked during the twenty five years 
since their inception. If instead of focusing on the rotation of the axes (here, as in Fig. O the thin straight 
lines) one focuses instead on the axis endpoints (small black dots), it immediately becomes apparent that 
the flat ribbon formed by the axes has, just like the parent Mobius strip, only one edge. Starting at the 
point labeled and proceeding counterclockwise along the edge (o ^ p ^ ... ^ v ^ w), after one 360° 
circuit around the C point (central large black dot) one reaches the point labeled 1. Continuing along a 
second circuit one returns to the point labeled 2, which, because it is coincident with the starting point 0, 
demonstrates that the axis ribbon has only a single edge. The single self intersection of the curve formed by 
the endpoints in (a), which is Fig. [2^ redrawn, suggests (but of course does not prove) the possible existence 
of a one-half-twist Mobius strip, whereas the three self intersections of this curve in (b), which is Fig. 
redrawn, suggests (but, again, does not prove - see Fig. [2b) the possible existence of a three- half- twist 
Mobius strip. 

Finally, we recall the well known maxim that the index of a path (here the circle ctq) is a property 
of the path. The index of the path becomes the index of the C point by a limiting process in 
which we shrink the path onto the point and observe that the index remains invariant under this 
transformation. L, satisfies this criterion also. 



FIG. 6: Index I^. (a) and (b) = +1. (c) Ij = —1. Here the gray ellipse facilitates following the 
retrograde rotation of the axes. The thick black markers and radial lines in (b) and in (c) relate to the line 
classification, and are discussed in Section III.D.2. 

What is the three dimensional structure generated by axis 7? As can be seen in Fig. [7] it is a 
segment of a cone, not a Mobius strip. This cone is analogous to the cones described previously 
that surround ordinary, i.e. nonsingular, points [14,15]. 




a 



b c 



a 



b 



c 



FIG. 7: 7 cones, (a), (b), and (c), respectively, are 3D views of the cones that produce projections (a), (b), 
and (c) in Fig. [6l 

B. Topological Indices of the Projection onto tq 

The classic construction of a Mobius strip, known to middle-school children everywhere, is as 
follows: take a strip of paper, twist one end through a half turn, and glue the two ends together. 
A similar construction serves to model the optical Mobius strips discussed here: attach short rods 
(representing ellipse axes a or (3) at their midpoints to a strip of flexible material, introduce one 
or three half-twists, and connect together the ends of the twisted strip. 

The strip can be twisted in one of two senses to form either a right- or left-handed screw. 
Rotating the ends of the strip to bring them together can also be done in one of two ways, when 
viewed from, say, above, either clockwise or counterclockwise. But the geometry and topology of 
the strip remain the same whichever choice (clockwise or counterclockwise) is made, and only the 
sense (right- or left-handed) and number (one or three) of the initial half-twists determines the 
nature of the strip. 

In the remainder of this report we drop for convenience the subscript "0" on tq, Xq, Yq, Zq, and 
Xo, its presence being understood. 

1. Indices Ta^p 

Following the (unfortunate) convention mandated for circularly polarized light, namely that the 
left-handed (right-handed) screw traced out in space by the rotation of the electric vector as the 
light propagates is labelled right-handed (left-handed) polarization, we attach a positive (negative) 
sign to left (right) handed Mobius strips. Adding in the number of half twists, we find four possible 
values for the twist index r, r = ±1/2, ±3/2. 

T can be obtained by following the rotation of the ellipse axes projected onto the rotating plane 
Tq. We label this projection T, and illustrate the procedure in Fig. [SI Here the angle C relative 
to the X-axis is measured as C ix) = arctan {Tz (x) > Tx (x))i where, Fig. [3ti, x measures position 
on (Jo. Cix) is unfolded (unwrapped) as required, and r is computed as r = A(/2tt, where 
AC = C(27r)-C(0). 

Although not an obvious geometrical or topological necessity, we find in Sq (and only in Sq) 
that in all cases the values of r obtained for axes a and (3 are always the same, even though the 
detailed geometry of the a and /3 Mobius strips may differ substantially. Accordingly, here we 
label this index Ta^p- In Section III.D.l we show that within the linear approximation only one- 
and three-half-twist Mobius strips are possible, so that generically r^^^ = ±1/2, ±3/2. 



FIG. 8: Index Ta,p. (a) Mobius strip with a single half-twist. This strip is the one shown in Fig. [T^. 
The origin of vector T is on ctq (gray circle) whereas its head terminates on the projections onto tq of the 
endpoints of the surrounding ellipse axes (small black dots) . T begins a 2tt circuit on cto at the the small 
white circle, and in the case shown here rotates clockwise through 180° to the large black dot that terminates 
the curve of endpoints. Here and throughout, this black dot is connected by a thin straight line to the initial 
white dot in order to emphasize that these two points terminate opposite ends of the same ellipse axis. The 
position of T relative to the X-axis is measured by angle C, = arctan (T^, Tx). At the start (end) of the 
circuit shown here C = 135° (C = -45°). The net rotation is therefore AC = -45° - 135° = -180°, and 
Ta.p = —1/2. (b) Mobius strip with three half-twists. This strip is the one shown in Fig. [T|d. Here T 
rotates counterclockwise through 4-540° over a 2tt circuit on (Tq, rotating by a little less than 360° over the 
lower loop and by a little more than 180° over the upper one. Ta^p is therefore Ta^p = -1-3/2. The starting 
point on ctq is arbitrary, and choosing a different starting point for the circuit changes the initial and final 
values of C, but not their difference A^. Although by convention we traverse ctq in the counterclockwise 
direction as viewed from the positive end of 7, the +z direction in Fig. [3^, traversing cto in the opposite, 
clockwise, direction does not change either the sign or magnitude of r^.p. This reflects the fact that the 
handedness of a screw — right or left — remains the same when the screw is viewed from any direction. 

a. Phase ratchet rules As T rotates its endpoint traces out one of the the curves shown 
in Fig. [HI The signed (plus for counterclockwise, minus for clockwise) crossings of this curve 
with the X and Z-axes can be used to determine by means of a minor variation of what we 
called the "phase ratchet rules" [31]. For a one-half-twist Mobius strip in a 2tt (47r) circuit on 
ctQ) Fig. Eh, the vector T rotates through vr {2tt), and the endpoint curve crossing sequence is 
XZ (XZXZ), or its cyclic permutation. For a three-half- twist strip T rotates through Svr (Gvr), 
Fig. [8)3, and the endpoint curve crossing sequence is XZXZXZ (XZXZXZXZXZXZ), or its 
cyclic permutation. An important complication, however is that the rotation of T is generally not 
monotonic, and T can, and often does, wander back and forth during its overall rotation, leading 
to additional crossings of the endpoint curve with the XZ-axes. Here the phase ratchet rules come 
into play: these rules state that adjacent terms of the same axis in the crossing sequence (these 
terms necessarily have opposite signs) are to be erased; thus, the sequence XZZZ for a 27r circuit, 
for example, is reduced to XZ, indicating that the Mobius strip has a single half-twist. These 
rules are illustrated in Fig. [9l 

2. Indices dTa^p 

As is evident from the phase ratchet rules, the complicated shape of the endpoint curve in 
Fig. Eh is not fully characterized by Tq,^^; here we discuss a second topological invariant that 
adds additional information about such curves. (A complete characterization of the curve would 
require, of course, an infinite set of indices or moments.) 




FIG. 9: Phase ratchet rules. (a) Endpoint curve. (b) Mobius strip. For the Mobius strip shown in 
(b) the endpoint curve in (a) is a comphcated figure with multiple crossings of the X and Z axes. As a 
result, the rotation of the vector T is so complicated that its net angle of rotation cannot be determined 
visually. The phase ratchet rules simplify the calculation and permit Ta^p to be obtained by inspection. In 
applying the rules one first lists the signed crossings of the endpoint curve with the X and Z axes in the 
order that they occur as one traverses the curve. Starting at the small white circle and proceeding along 
the curve, the sequence of signed crossings is Z+Z^X-X+X-Z-, where subscript + (— ) indicates that T 
rotates counterclockwise (clockwise) at the crossing. Applying the phase ratchet rules one erases adjacent 
terms in the same axis, in this case the terms Z+Z^ and X^X+, or to obtain X_Z_. From this 

canonical sequence it immediately follows, Fig. [Sk, that Ta^p ~ —1/2. The Mobius strip shown in (b) can 
be generated from an ordinary right-handed half-twist strip by grasping a segment of the ribbon between 
thumb and forefinger and twisting the segment by 180° (see the arrow in (b)). This introduces a pair 
of twists — one right-handed and the other left-handed. Although this operation changes the geometry of 
the strip, because the new right and left handed twists cancel the topology remains that of a right-handed 
one-half-twist Mobius strip. As discussed in Section IV, in a random optical field a surprisingly large 
fraction of all Mobius strips have complicated structures similar to the one shown here. 



A standard, widely used characterization of curves in the plane is the Poincare index. This 
index measures the winding number of the tangent to the curve, in our case the endpoint curve 
measured over one 27r circuit of co- Here we denote an analog of this index by dr, and cal- 
culate it as dr = arctan (dTz / ds , dTx / ds) using the finite difference approximation dr [xj) = 
aictan. {Tz {xj+i) — Tz {Xj) {xj+ij — Tx {Xj)) with uniform increment Asj = Sj+i — Sj = 
As > 0, where, Fig. [3^, s measures position on (Tq. 

Although the endpoint curves differ for axes a and /3, and although, like for t^^^, there is no 
obvious geometrical or topological necessity for dr to be the same for both axes, we find, just like 
foi" 'Tafi^ that in Sq (and again, only in Sq) in all cases dTa = dr^. Accordingly, in what follows 
we label this index dT^^p. 

Like Ta^f}, dTa^p takes on the values dTa^p = ±1/2, ±3/2. We show in Section IV that only ten 
of the sixteen possible combinations of r^,/? and dTa,f3 appear, so that these two indices are not 
completely independent one of the other. An example in which r^,/? = — 1/2 and dTa^p = ±3/2 is 
shown in Fig. [101 As shown in this figure, we take the sign of the local rotation of the tangent 
vector to be -|- (— ) if the vector rotates locally in the counterclockwise (clockwise) direction. 

But if dTa^p is even partially independent of Ta^fs, how can we be sure that dTa^p takes on only 
the half integer values ±1/2, ±3/2? The answer consists of two parts: in the first part, given 
below, we show that over a single 2tt circuit of o"o the vector tangent to the endpoint curve must 
rotate through nir, where n is an odd positive/negative integer. In the second part, which still 
remains to be accomplished, one shows that within the linear approximation |n| = 1,3. 

Over a 47r circuit of ctq, when we return to our starting point the tangent vector to the endpoint 
curve must also return to itself, so for a 47r circuit the tangent vector rotates through 2mr. But 



as illustrated in Fig. [TOl the endpoint curve for the 47r circuit has inversion symmetry, so that 
over any 27r circuit the vector rotates through rnr. The inversion symmetry of the endpoint curve 
reflects the fact that over a 47r circuit we visit the axis of each ellipse on ctq twice, arriving at 
opposite ends of the axis after a 27r circuit. From this inversion symmetry it follows that after 
any 27r circuit the tangent vectors at the beginning and end of the circuit must be antiparallel, so 
that n must be an odd positive/negative integer. 




FIG. 10: Tangent index dTa.p- As in Figs. [H] and H^, the central figure shows the curve (black dots) 
generated by the projection onto tq of the endpoints of axis a of the ellipses on the surrounding circle ctq- 
The point where ctq pierces tq is shown by the gray circle. The starting point of a 27r circuit on (Tq is shown 
by the white dot, the end of the circuit by the black dot that terminates the endpoint curve. The tangent 
vector to the endpoint curve is shown by arrows labelled with the orientation of the vector relative to the 
starting orientation, 0. As can be seen, the vector rotates through +37r over the circuit, so that dva^fs, 
which measures the net rotation of the vector, is dTa,f3 = +3/2. The inset in the upper left corner shows 
the endpoint curve for a 47r circuit. As can be seen, this curve is centrosymmetric about ctq- ^0,13 can be 
determined by inspection using the phase ratchet rules. Here the sequence of crossings of the curve with the 
XZ-ax.es is Z^X^Z^Z-^^; after applying the rules this reduces to showing that Tq.^ — —1/2. 



3. Index dr^ 

As shown in Fig. [71 axis 7 does not generate a Mobius strip. It may therefore appear surprising 
that as Tq rotates along cto the projections onto this plane of the 'y-axis endpoints generate a two 
turn helix. This helix, which can be either right- or left-handed, is characterized by tangential 
winding number dr^, the equivalent for axis 7 of the index dr^^i^ for axes a and /3. As can be 
seen in Fig. \7\ the endpoints of axis 7 return to themselves after one 27r circuit, so dr^ = itn; 
we find in all cases \n\ = 2, reflecting the fact that both turns of the helix are identical. Like a 
Mobius strip, a helix retains its handedness, and therefore its winding number, when viewed from 
any direction. Examples of dr-y are shown in Fig. [TT\ which also shows that the index T~y, the 7 
equivalent of Ta,/3, is always zero, and is therefore of no interest. 

C. Topological Indices of the Projection onto ttq: Indices tTq jj, dna^p, and dTr-y 

The projections of the axes of the ellipses on (Tq onto the plane ttq in Fig. [3b generate three 
winding numbers that are analogs of T^^p, dTa^p, and dxy. 



FIG. 11: Tangent index dr^. Axis 7 endpoints (black dots) are plotted as a function of arclength s = rox, 
where x measures position on co, Fig- and is the radius of (Jq, (a) The endpoints generate a left- 
handed helix. Starting at the white dot, the tangent vector to the curve (arrow) rotates monotonically in 
the positive, counterclockwise direction through 47r, and dr^ = +2. (b) Here the endpoints form a right- 
handed helix, the tangent vector rotates through 47r in the negative, clockwise direction, and dr-^ = —2. (c) 
Index Tj. Vector T oscillates back and forth with a net winding angle of zero. Starting at the white dot 
and proceeding clockwise, the crossing sequence of the endpoint curve with the XZ-axes is Z-Z^Z-Z^. 
Applying the rules of the phase ratchet yields an empty sequence, so T-y = 0; this is the value of for all 7 
cones. Unlike a and /3, whose endpoint curves are typically complex. Figs. [9l and fTOl 7 endpoint helices 
always have smooth, simple shapes. 

The ttq analog of Tq,^^ we label ira,f3', like Ta^f3, i^a^p takes on the values ±1/2, ±3/2, is the same 
for axes cx and /3 in Sq, and can be determined by inspection using the phase ratchet rules with 
axes XZ replaced by Y Z. All 16 combinations of r^^p and TTa,p appear in our simulations, showing 
that these two indices are independent one of the other. 

The vTo analog of dT^^p we label diTa^p', like dr^^fs, d-Ka,p takes on the values ±1/2, ±3/2, is the 
same for axes a and /3 in Sq, and can be determined by following the rotation of the tangent vector 
to the endpoint curve in vro. We show in Section IV that, just like for r^,/? and dTa^p, only ten 
of the sixteen possible combinations of TTa,^ and dira^/s appear, so that these two indices are not 
completely independent one of the other. 

The ttq analog of dr^ we label dir^; like dr^, dn^ takes on the values ±2, and can be determined 
by following the rotation of the tangent vector to the 7 endpoint curve. All 4 combinations of 
dr^ and dny appear in our simulations, showing that these two indices are independent one of the 
other. 

7ro,/3, d-Ka,^-, and dvr^ are illustrated in Fig. [TH 

D. Line Classification 

In addition to determining winding numbers la^p for axes a and /3 (Section III.A.l), and for 
axis 7 (Section III. A. 2), the projections of these axes onto Eq determines another important prop- 
erty of ellipse fields: the line classification. In this classification, introduced by Berry and Hannay 
[18], one counts the number of straight streamlines that terminate (or equivalently originate) on a 
singularity; here we denote this number by A. We discuss this classification for C points for axes 
OL and /3, and we then extend it to axis 7. 



FIG. 12: Indices tTq^^, diTa^p, and diTj. (a) TTa^p and d'Ka.fi- '^a.f3 measures the net rotation of vector P. 
Starting at the small white circle, the crossing sequence of the endpoint curve (small black dots) with the 
FZ-axes is Y^Z^Z+Z^, which, after application of the phase ratchet rules reduces to Y^Z^, showing that 
TTa,/? = —1/2. diTa.p measurcs the net rotation of the tangent vector (small arrow) to the endpoint curve. 
This vector is easily seen to rotate by Stt in the negative, clockwise direction, and dna.p = —3/2. (b) dir-y. 
d-Kj measures the net rotation of the tangent vector (small arrow) to the 7 helix, plotted here, as in Fig. 
11, vs. arclength s on ctq. As can be seen, the tangent vector rotates in the negative, clockwise direction 
by 47r, and dir^ = —2. The data in (a) (in (b)) for 7Ta,f3 and d'Ka.p (for c?7r^) are for the same Mobius strip 
(same 7 cone) as that shown in Fig. [TO] (in Fig. [TTk). 

We start be reviewing well known results about application of the line classification to C points. 
Berry and Hannay showed that within the linear approximation, for umbilic points there are only 
two possibilities A = 1 and A = 3 [18]. At an umbilic point the eigenvalues of the curvature tensor 
(matrix) become degenerate. This tensor can be represented by an ellipse, and at an umbilic 
point the major and minor axes of the ellipse become equal and the ellipse degenerates to a circle. 
From this it easily follows that the results of Berry and Hannay are applicable to C points, which 
as discussed in Section II. A, are points where the eigenvalues of the coherency matrix become 
degenerate and the polarization ellipse degenerates into a circle — the C circle. 

When the path that encircles the C point is a circle with sufficiently small radius, A is the 
number of lines in plots such as those in Figs. [2l HI and [6] that are radially directed towards the C 
point. Labelling A by the relevant axis we have from Berry and Hannay, Aq, = A^ = Aq,^^ = 1,3. 

The reason why Aq, = A^ for a C point is as follows: An important geometric property of 
the lines in Figs. 2, 4, and 6, is that lines on opposite sides of gq are orthogonal. Thus, for 
every line that points towards the C point there is a line on the opposite side of the circle that is 
tangent to (Tq. But as discussed in Section II. A. 1, for o"o sufficiently small the lines representing 
the projections onto Sq of axes a and /3 are also orthogonal, so that wherever a is tangent to cjo, 
/3 points towards the C point, and vice versa. 

Here we determine A^^^ numerically, finding the number of points on do at which the orientation 
9 of the ellipse axis, measured counterclockwise relative to the x-axis, equals x, the angular position 
of the axis center on do. Fig. [3l Specifically, calculating 6 for say axis ct as = arctan (a^, a^), 
where {ay) are the components (projections) of a along the x- and y-axes, respectively, we solve 
^ (x) = X + '^'''"1 where n = —4 ... +4, as required [32, p. 258]. Fig. [TBb illustrates a graphical 
implementation of this method. We note that an alternative, analytical, method would be to use 
the discriminant given by Dennis [33], writing the Stokes parameters in Sq in terms of ax and ay 
[13], but the resulting expressions for the general case are so unwieldy as to be impractical. 

There is an important connection between I^^p and Aq,^^. Because la^p = +1/2,-1/2 and 
A-a,/3 = 1,3, it might be thought that there are four possible combinations of these indices. Berry 



and Hannay showed, however, that this is not the case and that the number of possible combinations 
is only three: (i) la^p = +1/2, ^cx,p = 1) the lemon] {ii) I^^p = +1/2, ^a,p = 3, the monstar; and 
(iii) lafi = —1/2, ^a,p = 3, the star. These three combinations are illustrated in Fig. [T3]b — f. 
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FIG. 13: Aq_^. (a) Graphic implementation of the numerical method described in the text for determining 
A. The solid curve is the angle 9 (x) that the projections of axes a make with the x-axis. A counts the 
intersections (small white circles) of this curve with the dotted lines 9 (x) = x + nir, where here n = —1, —2, 
and —3. These data are for the projection shown here in (c). (b) A Itq^^I — 1/2 lemon. Here la.jS — +1, 
Aa,/3 — 1, and Tq,^ = +1/2 . Throughout, a thick line marks an axis that points directly at the C point 
(central dot); perpendicular to this is an axis, marked by a white line, that is tangent to the cto circle, (c) 
A |Ta,/3| = 1/2 star. Here la.fi — —1, A^^.p — 3, and t^^/j = —1/2. (d) A |tq_;3| = 3/2 star. Here la^p = — 1, 
Aa,/3 = 3, and Ta,(3 = +3/2. (e) A \Ta.j3\ = 1/2 monstar. Here Ia,p — +1, Aq,^^ = 3, and Ta^p = —1/2. (f) 
A |tq./3| = 3/2 monstar. Here Ia,i3 — +1, Aa,p — 3, and Ta.p = +3/2. As discussed in the text, there are 
no \Ta.j3 \ = 3/2 lemons. 

There are also important connections between Aq,^^, and Ta^/s, i^a^p, the phase ratchet rules, and 
the endpoint curves in Figs. [SI and [TUl these connections follow from the facts that when the 
endpoint curve of say cx crosses the Z-axis of the plane tq, the projection of a onto Eq is tangent 
to (To, whereas when the endpoint curve crosses the Z-axis of ttq, the projection of a. onto Sq is 
radial (along the radius of do). These important geometrical results are illustrated in Fig. [TH 

We therefore have the following: 

(i) Endpoint crossings with the Z-axis in tq and the Z-axis in vro occur at opposing points on 
(To (points separated by 180°). 

(ii) ^a,f) equals the number of endpoint crossings of the Z-axis in tq: at each crossing the axis 
projection onto Sq is tangent to (Tq, implying the presence of a radial projection at the opposing 
point on itq; by definition, A^^p is the number of such radial projections. 

{iii) Aa,/3 also equals the number of endpoint crossings of the Z-axis in ttq, because at each such 
crossing the axis projection onto Sq is radial. 
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FIG. 14: Z-axis crossings. Ellipse axis centers (endpoints) are shown by small gray (white) circles, whereas 
the axis itself is shown by the connecting thick line. The projected endpoint curves in the planes tq (the 
X Z-p\ane) and ttq (the FZ-plane) are shown by the connected black dots, (a) When the endpoint curve in 
To crosses the Z-axis, the ellipse axis lies in the yZ-plane, and its projection onto Eg is tangent to the erg 
circle, (b) When the endpoint curve in ttq crosses the Z-axis, the ellipse axis lies in the XZ-plane, and its 
projection onto Sq is radial. 

From (ii) and (iii) follows that the number of endpoint crossings with the Z-axis is the same 
for both To and ttq, and that this number must be either 1 or 3. Below, for simplicity we refer to 
an endpoint crossing with the Z-axis as a "crossing" , and we have the following: 

(/) If Aa^f3 = 1, \Ta,f3\ = 1/2: from the phase ratchet rules if |tq,^^| = 3/2 there must be more 
than one crossing. Thus, all lemons are Mobius strips with a single half-twist. The converse, 
however, is not true, and a one-half-twist strip can have Aq,^^ = 3, and can therefore be a star or 
a monstar. Figs [TSb.d 

(//) If |tq,^^| = 3/2, then Aq,^^ = 3, so all Mobius strips with three half-twists are either stars 
or monstars. This important result is to be credited to Dr. Mark R. Dennis (Bristol), who, using 
the methods of differential geometry, was first to derive it [34]. 

(///) Within the linear approximation the Mobius strips surrounding C lines can have one or 
three half-twists — strips with say five half-twists cannot occur within this approximation. Strips 
not necessarily associated with C lines, or with L lines (which have four half- twists), are also 
possible; an example is the two-half-twist (one- full- twist) Mobius strips that surround ordinary 
ellipses [14, 15]. 

2. Aj 

For the line classification for axis 7 we find the following: If = +1, A^ = 0,4, whereas if 
= —1, A^ = 4; examples of all three cases are shown in Fig. El where, as in Fig. [13] for A^^^, 
axes that point toward the C point are marked by thick lines. Unlike the case of axes a and 
/3, however, at opposing points on ctq, 7 axes are not perpendicular, but are parallel. Whether 
A^ = 0, or A^ = 4, depends on whether the endpoint curve in ttq crosses the Z-axis, because the 
geometrical considerations shown in Fig. 114b hold for all three axes — a, /3, and 7. In Fig. [T^] 
we show the three endpoint curves that correspond to the three cases in Fig. [6l As can be seen, 
when A^ = ( A^ = 4) the endpoint curve crosses the Z-axis (4) times. 
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FIG. 15: Axis 7 endpoint curves in ttq. Panels (a), (b), and (c), are the endpoint curves corresponding to 
panels (a), (b), and (c) in Fig. El respectively. As can be seen, here in (a) the endpoint curve does not 
cross the Z-axis and in (a) of Fig. [6] A^ = , as expected. In (b) and (c) there are four crossings and in 
(b) and (c) of Fig. =4, again as expected. 



E. Index Summary 

We summarize here the rather complex results discussed above for the various indices. All 
indices are the same for axes o: and (3 and are therefore labelled "a,/3". 

(z) Index /^^^ measures the rotation of the projections of axes a and (3 onto Sq, the plane of 
the C circle, Fig. [2l Generically /q^,/3 = ±1/2. A positive (negative) value for this index, and 
for index below, implies that the axis projections rotate in the same (opposite or retrograde) 
direction as the path on which they are measured. 

(ii) There are either one or three points on a sufficiently small circular path centered on the C 
point at which the axis projections point directly towards the point. The number of such points 
is measured by the line classification index A^^^. Unlike Ia,f3y ^a,i3 is not a topological invariant, 
but is nonetheless an important, highly useful characterization. 

{in) If Ia^f3 = +1/2 and Aa,i3 = 1 i^a,p = 3), the C point is dubbed a "lemon" (a "monstar"). 
Fig. [131 Lemons always correspond to Mobius strips with a single half- twist, whereas monstars 
can be Mobius strips with one or three half-twists. If la^/s = —1/2, then A^^^ = 3, and the C 
point is dubbed a "star". Stars can correspond to Mobius strips with one or three half-twists. 
Conversely, if a Mobius strip has a single half-twist without any mutually cancelling additional 
twists (Fig. O)) it must be a lemon, whereas if it has three half-twists it can be either a star or a 
monstar. The above relationships are the same for both right- and left-handed Mobius strips. 

(iv) Index 1^ measures the rotation of the projection of axis 7 onto Sq; generically I-y = ±1, 
Fig. m The line classification index A^ measures the number of points on the path at which 7 
axis projections point directly towards the C point. If Xy = +1, A..y = 0,4, whereas for = — 1, 
A^ = 4. Ia,/3 and I-y are independent, and all four combinations of these indices can appear. I-y 
and Ay are independent of the number of twists or handedness of the Mobius strip. 

(v) The magnitude of index Ta^p measures the net number of half twists of the Mobius strip, 
the sign of Ta^fj measures the handedness of the strip (+ = left-handed, — = right-handed). Figs. 
Elandlll There are four possible values for r^,/?: = +1/2, —1/2, +3/2, —3/2. 

(vi) There are in principle 2 x 2 x 4 = 16 possible triple combinations of /q.,/3, A^^^, and r^^p, but 
(see {in) above) only 10 are allowed: la^p = +1/2, A^^^ = 1, |tq^/3| = 1/2 (a one-half-twist lemon); 
Ia,/3 = +1/2, Aq^^ = 3, \Ta,/3\ = 1/2,3/2 (a one- and a three-half-twist monstar); = —1/2, 
Ao,/3 = 3, Itq^^I = 1/2,3/2 (a one- and a three-half-twist star). 

(vii) Ta^p is independent of and A^, and for each of the three allowed and A^ combinations 
in [iv) above, all four values of Ta^p appear, so that in total there are 12 triple combinations of 1^, 



A^, and T^^f^- 

(viii) Whereas Tq,/? measures the net signed number of times axes a and f3 rotate or loop 
around the Une of the path, index measures additional signed rotations or loops that do not 
encircle this line, Fig. [lOl dTa,/3 takes on the four values dra^jj = ±1/2, ±3/2. A right-handed 
one-half- twist Mobius strip (tq,,/? = —1/2), for example, can have dTa^p = ±3/2 due to additional 
loops that are "up in the air" , as it were. 

{ix) Indices TTa^p and d'7Ta,/3 measure oscillations of axes a and P that form signed loops not 
captured by indices Ta^p and dTa^p'- T^a,f3 = il/2,±3/2; diTa^p = ±1/2, ±3/2, Fig. [T21 7ra,i3 is 
independent of diTa^/s is independent of dTa,p- 

(x) Indices dr^ and (iyr^ measure "up in the air" oscillations of axis 7 that form signed loops 
which do not enclose the line of the path. Fig. [TTJ These indices describe the cones generated by 
axis 7, Fig. [71 are independent of one another, and of the corresponding Q;,/3 indices, and take the 
values: dr^ = ±2; dvr^ = ±2. 

Altogether there are a total of ten indices that characterize C points, their Mobius strips, and 
their cones: Ia,p, I'y, -^0,13, ■^'y,Ta,i3, dTa,/3,TTa,/3, diTa,/3,dT^,d7r^. If independent, these indices could 
generate 2^ x 4"^ x 2-^ = 16384 different C points (C lines). The connections discussed above reduce 
this number to 7680. In Section IV we discuss additional connections that reduce this number 
even further, and present statistical data for the various index combinations. 

IV. STATISTICS 

Here we list the most important statistical properties of the geometrically and topologically 
distinct Mobius strips that appear in Sq. We emphasize that these statistics change importantly 
when the plane of observation S is rotated away from Sq by small, but finite angles. 

As noted in Section III.E, the ten indices Ia,i3, 1-y, ■^a,i3T^-yiTa,f3i dTa^i3,Tra,f3, diTa^f}, dr^jdir^, if 
independent, could generate 16384 different C points (C lines). This number is reduced to 7680 
by the constraints discussed in Section III.D. In this section we search for additional constraints 
(selection rules) that further reduce the number of allowed index combinations, and we present 
statistical data for the probabilities of the remaining combinations. 

We search a database containing 10^ independent, randomly generated realization, after remov- 
ing all realizations in which any one of the ten indices differed by more than ±0.001 from an integer 
or half integer. This eliminated ~ 2% of the realizations. The the remaining index combinations, 
which all obey the constraints discussed in Section III.D, have a very broad distribution. This is 
illustrated in Fig. [16] using a modified Zipf plot. As can be seen from this plot, it is unlikely that 
the number of configurations appearing in an arbitrarily large database would substantially exceed 
1,200. 

A. Selection Rules 

It is useful to formulate constraints among indices in the form of selection rules — rules that 
state which combinations of indices are forbidden. Here we list 13 such rules: 8 rules involve pairs 
of indices — binary rules; 5 rules involve index triplets — ternary rules. 

The search for binary (ternary) rules is facilitated by 2D (3D) plots such as those shown in Fig. 
[T71 There are, in principal, a total of 45 2D plots and 120 3D plots that need to be examined 
and interpreted. Our algorithms, however, display for analysis only those plots in which there are 
missing indices, reducing the number of plots to 8 for 2D, and 62 for 3D. Of the 3D plots, 57 do 
not yield new ternary rules, because the absence of configurations results from combinations of the 
binary rules. An example of this is shown in Fig. [T7b . 




FIG. 16: Modified Zipf plot. The logarithm of the number of exemplars, N^ank, found for a given ten- 
index configuration is plotted vs. the rank of the configuration, where configurations are ranked in order 
of decreasing probability. The straight line is N^ank — Noexp {—K x Rank), with No — 4950 and 
K = 0.00575. The tabular inset lists the four most probable configurations (the four unresolved points in 
the upper left corner). There are some 2 x 10** exemplars of each of these four configurations in the complete 
10^ entry database. As can be seen, these four most probable configurations are all three-half-twist stars. 



Altogether, the 13 selection rules reduce the number of index configurations from 16, 384 to 
2, 104. This is almost twice the number observed. Fig. \T6\ indicating the possible presence of 
additional, higher-order rules. There is no simple way to systematically search for such rules, 
however, since in the overwhelming majority of cases the absence of index configurations for, say 
the 210 quaternary combinations, will be due to binary and ternary selection rules acting together. 
Because of this we have not attempted higher-order searches. 




FIG. 17: Index combinations, (a) 2D plot for binary index combinations. (b),(c) 3D plots for ternary 
combinations. In each case empty grid points correspond to missing, and therefore presumed forbidden, 
combinations. The areas of the black circles are approximately proportional to the relative probabilities of 
the different combinations, (a) Binary selection Rule 2.5 is derived from this plot, (b) Plot from which 
ternary selection Rule 3.4 is derived, (c) Here all missing indices are due to combinations of binary rules. 
The interested reader may enjoy analyzing the 16 configurations (marked by x's) missing from this plot in 
terms of the rules. 



1. Binary Rules 



Listed below are the 8 binary selection rules (forbidden index combinations). The number 
of configurations each rule eliminates when acting alone is listed in parentheses. The number 
of configurations eliminated by combinations of rules may be difficult to tally because the rules 
interact in complex ways. Together, these 8 binary rules reduce the number of index configurations 
from 16, 384 to 2, 712. 



Rule 2.1. Ia,i3 
Rule 2.2. I ' 
Rule 2.3. 
Rule 2.4 
Rule 2.5 
Rule 2.6 
Rule 2.7 
Rule 2.8 



7 ~ 



-1/2 & K^p = 1.(4096) 
-1 & = o'. (4096). 
= 3/2 & A„,^ = 1. (4096) 
= 3/2 & Ao'/3 = 1. (4096) 
= 3/2 & dT^^p / T^^p. (6144) 
= 3/2 & d7r„,^ / T:o,,p. (6144). 
= 3/2 & sign(dT^) ^ sign(r«,^). (4096) 
= 3/2 & sign((i7r^) ^ sign(7r„,/3). (4096). 



2. Ternary Rules 

The five forbidden ternary selection rules listed below acting together (without the binary rules) 
reduce the number of index configurations from 16, 384 to 8, 064. 

Rule 3.1. Ia,f3 = 1/2 & \Ta,p\ = 3/2 & Ivr^,,/?! = 3/2. (2048) 
Rule 3.2. Ia,f3 = 1/2 & A^ = & |t„,^| = 3/2. (2048) 
Rule 3.3. iJi, = -1/2 & A^ = & |r„,/3| = 1/2. (2048) 
Rule 3.4. = -1 k |r^,^| = 3/2 k = 3/2. (2048) 

Rule 3.5. A^ = & sign'((ir-y) / sign(T„'/3)- (2048) 

These five rules cover all configurations not observed in the 62 3D plots discussed above. 

Except for Rules 2.1—2.4, which are based on known C point results or easily verified geometrical 
constraints, the above rules are provisional, and await confirmation from analysis. These rules, 
however, are practical rules, since any configuration not seen in 10^ samples in unlikely to appear 
in an experiment. The fact that the rules reduce the number of configurations only to 2104, which 
is nearly twice the number actually observed. Fig. [TBI may imply the existence of additional rules 
or may be due to other, presently unknown, factors. 



B. Probabilities 

Listed below are the probabilities of occurrence of various configurations in our database. We 
emphasize that these are not densities, and therefore do not directly correspond to the quantities 
that would be measured in an experiment. In order to compute densities each realization needs 
to be weighted by a Jacobian (with units of area~^) that, like the realization itself, is a function 
of the wave field parameters in Eq. ([3j) . 

At present there are two Jacobians available [9]. Defining J = V (Re {E • E)) x V (Im [E • E)) 
these are 

J2D = \ Jz\ , (7a) 
JsD = \J\. (7b) 



In Eq. (j7ap the z-axis is normal to the plane of observation, here S. J2D is suitable for counting 
2D densities in the plane, such as the density of C points that appear where C lines pierce the 
plane. J^d measures 3D quantities, such as the length of C line/unit volume. But neither of 
these two Jacobians, at least as written, depends sensitively on the angle between the S and Sq. 
This sensitivity is essential, because as emphasized above, the statistics change importantly when 
E rotates away from Sq by a small, but finite angle. 

The probabilities listed below, although not directly connected to experiment, are, in principle, 
calculable from theory. We note that eliminating the Jacobian from the theory, in fact, greatly 
simplifies the calculation. 

Among the endless probabilities P that can be extracted from our database we believe the 
following to be of greatest interest: 

(1) Positive and negative values for T^^p and vr^^^ occur with equal probability P, 
and P (t„,;3 = ±\) =P (7r„,;3 = ±\) = 0.80, P {jg,f} = if) = P (vr„,/3 = ±1) = 0.20. Through- 
out, we use the notation r^^p = zt^, tTq^^ = ±i, instead of |t„^^| = |7ra,/3| = 1/2, etc., to emphasize 
that in searching our database we search for each sign combination separately. When in later 
tables the results are sign dependent we list each sign combination. 

(2) The frequency of occurrence of Mobius strips with one and three half-twists that are lemons, 
monstars, and stars, are listed in Table 1. The selection rules forbid lemons with r^^p and T^a,[) 
equal to ib|. 



Table 1. Lemon, monstar, star, and indices Tq,^^ and T^afi- 





Lemon (I^^ = i,A^^ = 1) 


Monstar (/^^ = ^, A^^ = 3) 


Star(/^., = -i,A„, =3) 


^\ 


0.39 


0.11 


0.30 


±1 





0.005 


0.20 



(3) Combinations of Ta,f5 and i^a^p are listed in Table 2. As can be seen, it is some four times 
more likely that these two indices are the same than that they differ. 

Table 2. Indices r^,/? and tTq^^. 



Indices 


= ±^ 






0.17 


0.026 


^2 


0.026 


0.025 



(4) Table 3 lists probabilities for combinations of Ta^p and dTa^p. These probabilities are the 
same for tTq^/j and dira^p combinations, and are sign dependent. Combinations with zero probability 
are forbidden by selection rules. For the case |tq^^| (|7rQ^/3|) = ^, the number of occurrences in 
which dTa^/3 {d-Kafi) 7^ Ta^p i'^a,/3) excccds those in which (ir„^^ {d-Ka^f}) = (tTo^/j), reflecting the 
fact that somewhat more than half of all Mobius strips have a complex structure containing an 
additional, self cancelling pair of half-twists, (see Fig. [9]). 



Table 3. Indices Ta^/s (vTo,/?) and dra^jj {dna^/s). 



Indices 


(^Q./S) = +h 


1 

2 


^2 


3 
2 


dTa,i3 {dTTa,!}) = +\ 


0.19 


0.11 








1 
2 


0.11 


0.19 










0.09 


0.008 


0.10 





3 
2 


0.008 


0.09 





0.10 



(5) Correlations between the orientations of axis 7 and axes a, (3 are reflected in Tables 4 and 
5. Worthy of note is that the possibly surprising statistical equivalence of Tq,^^, and tt^^^ evident 
in all the other tables is broken in Table 5 for A-y. 



Table 4. Indices and r,,.,i (ttq,^/?). 



Indices 


Ta,l3 (TTa,/?) = ±^ 




I^ = +l 


0.16 


0.086 


-1 


0.23 


0.014 



Table 5. Indices Aq,^^, A^, Ta,/3, and tTq^/j. 



Indices 




^2 




^2 


Aa,/3 = 1 


0.19 





0.19 





3 


0.21 


0.10 


0.21 


0.10 


A^ = 


0.073 


0.073 


0.10 


0.042 


4 


0.33 


0.028 


0.30 


0.059 



V. SUMMARY 

Prior studies of C points and C lines have concentrated on a single projection of the major and 
minor axes of the ellipses surrounding a C line onto a plane, in most instances the plane we call Sq, 
the plane of the C circle. Examining the full 3D arrangement of the surrounding ellipses we found 
that their major and minor axes, ot and /3, generate Mobius strips, whereas the normals to these 
ellipses, y, generate a section of a cone. The Mobius strips have either one or three half-twists, 
and form segments of either right-handed or left-handed screws. 

The So projections of a and f3 give rise to two well-known indices [1,3]: the conserved topo- 
logical index 7q,/3) and the highly useful, albeit nontopological, line characterization A^,/?. These, 
and all other indices in Sq, have the same values for axes a and (3, whence the subscript a,/3. 

Ja.g measures the rotation of the projections of a, j3 around the C point, and takes on the generic 
values la^p = ±1/2. On a circular path centered on the C point, here labelled ctq, A^^^^ counts 
the number of times an axis projection points directly at the C point. Generically, A^^^j = 1,3, 
and there are three possible combinations of /q,^^ and A^ ^g: the lemon (/a,/3 = 1/2, A^^^ = 1); the 
monstar (/q,/3 = 1/2, A^,^^ = 3); and the star {Ia,i3 = ~l/2, A^^^ = 3). 

We introduced 8 new indices. 

Two of these involve the projections of axis 7 onto Sq. As is the case for axes a. and /3, for 
axis 7 this projection gives rise to a conserved topological index, I^, and a line characterization, 
A^. Generically = ±1, and A^ = 0,4, and also here there are only three possible combinations: 

= 1, A^ = 0, 4; = -1, A^ = 4. 

Six new indices arise from two new projections: the tq projection and the ttq projection. 

The major and minor axes of the surrounding ellipses generate two winding numbers for each 
projection: r^^p and dTa^p for the tq projection, iTaS a-^d '^'^a,/3 for the ttq projection. Each of these 
four indices can take on four different values: ±1/2, ±3/2. T^^p measures the number of half-twists 
of the Mobius strip, the other three indices measure other, more subtle structural features. Ta,j3 
and T^a,p a-i's independent, and all 16 combinations of these indices are found. Tq,^^ and dr^^/s, and 
also and dua^p-, are not completely independent, and only 10 of the 16 possible combinations 
of each index pair are allowed. 



The To and ttq projections of axis 7 generate the last two of our 8 new indices. These are dr-y 
and dvr^, each of which can be ±2. 

A set of simple rules, the phase ratchet rules, were formulated that permit Ta^p and T^afi to be 
determined by inspection from the tq and ttq projections; dTa,(} and diTa^f} can also be determined 
by inspection from these projections. 

Of the 16, 384 combinations of indices that would be possible in the absence of all restrictions, 
some 1, 150 were observed in a database containing 10^ randomly generated realizations. This 
number of distinct C lines exceeds by a factor of ~ 350 the three types of lines (lemon, monstar, 
star) recognized previously. Thirteen selection rules were formulated that reduce the number 
of possibilities to 2, 104. These rules include the well-known restriction that gives rise to the 
lemon-monstar-star trio, together with 12 new rules; the most important of these new rules are: all 
lemons are one-half-twist Mobius strips, and [34], all three-half-twist Mobius strips are either stars 
or monstars. The converse of the lemon rule, however, does not hold, and one-half-twist Mobius 
strips can also be monstars or stars. 

Statistics of the most important configurations were also presented. Most notable of these is 
that 80% of all Mobius strips have a single half-twist, the remaining 20% three half-twists. Also 
noteworthy is the fact that although all three-half-twist strips are monstars or stars, the majority 
of monstars and stars are, in fact, one-half- twist Mobius strips. 

A recurrent theme was the complexity of the various Mobius strips; this complexity was illus- 
trated by numerous examples, and described geometrically. Because of this complexity, a more 
complete mathematical description of these objects, which includes the dependence of the various 
indices, and their selection rules, on the parameters of the optical field, may be rather complicated. 

Recent advances make experimental measurements of the Mobius strips and cones feasible in 
both the microwave and optical regions of the spectrum. Such experiments would do much to 
further our knowledge of the structure of real, as opposed to simulated, 3D optical fields. 
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